Ultra-short solitons and kinetic effects in nonlinear metamaterials 
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We present a stability analysis of a modified nonlinear Schrodinger equation describing the propa- 
gation of ultra-short pulses in negative refractive index media. Moreover, using methods of quantum 
statistics, we derive a kinetic equation for the pulses, making it possible to analyze and describe 
partial coherence in metamaterials. It is shown that a novel short pulse soliton, which is found 
analytically, can propagate in the medium. 

PACS numbers: 42.25.Bs, 42.25.Gy, 78.20.Ci 

As early as in 1945, Mandelstam studied the electromagnetic properties of materials exhibiting a negative index of 
refraction P, Q . The concepts of negative permittivity e and permeability [i and the consequences of negative group 
velocity were moreover given attention by Pafomov jj] and Agranovich and Ginzburg Q , and have since then been 
discussed by several authors, e.g. 0, IS 0] • Recently, such materials have been produced [1JE3, and experimental 
verifications of previous theoretical issues are currently attracting much interest (see, e.g. |lll Il2j | for a review). 
From a sharp resonance in the material response to the applied external field, one may obtain negative e and /i. The 
normal procedure for obtaining negative-index-of-refraction materials is to put together two structured materials that 
both have negati ve p ermittivity and negative permeability, such that the resulting composite material has a negative 
refractive index Il2l . Since such materials may have a nonlinear response as well, a nonlinear metamaterial can be 
constructed [l3l.ll4j 

In the present Brief Report, we shall consider ultra-short electromagnetic pulses in left-handed materials |l4|j. We 
investigate the stability properties of a pulse, and show that the metamaterial can influence its dynamics in significant 
ways. Moreover, the effects of partial pulse coherence is studied using the Wigner function approach. It is also shown 
that the dynamical equation for modulated short pulses admits a novel localized envelope soliton solution, which 
should be experimentally verified. 

Given the electric field E(z, t) exjp(ikoz — iuj^t) of a short electromagnetic pulse, the evolution of the pulse envelope 
E in the slowly varying envelope limit, i.e. ko ^> (d z — 2(<9fc /<9wo)<9 t ), is governed by (e.g. [Tip 
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where we have introduced the parameters 

X^[n/2V g n 2 — (7 + (i)/2n]. Moreover, f3 — 2ncj a /uj p , uj p is the plasma frequency, x ta) is tne third order suscep- 
tibility of the medium, —a is the group velocity dispersion (GVD), V g — (dko/dujo)" 1 is the group velocity, fi is the 
magnetic permeability, n is the refractive index, 7 = d(u>on) / 'du>o, and the intensity is given by / = \E\ 2 . 
There is a stationary solution E(t, z) — Eq exp(iAfcz) of Eq. lfT|l. corresponding to a nonlinear phase shift 



Afc = «i(l - k 2 I )Io, 



(2) 



where Iq = E 2 . Linearizing around this stationary solution, according to E(t, z) = [Eq + Exit, z)\ exp(iAfcz), where 
I .Ei I < E , Eq. (TJ gives 
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By making the ansatz E\ — E\ r + iEu, and assuming a harmonic dependence in the form E± r , En oc exp(iKz — iClt), 
we can split Eq. @ into its real and imaginary parts, and obtain the dispersion relation 



K = -2ani ± n 
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1/2 



(4) 



2 



Letting K — — 2erf2io — iT, where T > 0, we thus have the modulational instability growth rate 



r = n 



-,1/2 



(5) 



The effects of the material parameters a, k±, and a on the modulational instability growth rate 10 can be divided 
into a number of cases. We note that in the coherent case the effect of the derivative nonlinearity is always to suppress 
the modulational instability. 

First, in the case of positive GVD, i.e. a < 0, there is no growth in a focusing medium (x^ > 0; n± > 0) if the 
higher order nonlinearities through K2 are neglected. If the quintic nonlinearities are taken into account for focusing 
media, we see from JSJl that a positive growth rate is possible, given that k±K2 > 0. Since K%/k 2 = 2f3n, we will thus 
have growth for normal and focusing media, where both the third order susceptibility and the refractive index are 
positive. In focusing metamaterials with positive GVD this is no longer possible, and the modulational instability thus 
disappears. If we instead have a defocusing medium, i.e. x^ < 0, this implies K\ < 0, and an\ > 0. Thus, growth 
is possible even if K2 — 0. For normal bulk media, <xk\K2 < which increases the modulational instability growth 
rate . However, for left-handed metamaterials we have uk,\K2 > 0, which implies a suppression of the modulational 
instability. 

Secondly, there is the case of negative GVD, so that a > 0. If we have a focusing medium, so that x^ > and 
K\ > 0, the effect of a metamaterial in which K2 < 0, is to increase the growth rate J5J, while normal media tend to 
inhibit the modulational instability. However, if the medium is defocusing, the coherent modulational instability is 
not possible in a normal medium, while in left-handed metamaterials we can obtain a positive growth rate from (jSJ) . 

Thus, we see from the cases above that the interplay between the material coefficients is complex, and requires 
careful consideration when studying the stability properties of waves in metamaterials. 

In order to take the effects of partial incoherence into account, we define the space-time correlation function for the 
electric field as F(z + , z_, t+, t-) = (E*(z+,t+)E(z-,t-)) > where z± = z±£/2, t± = t± r/2, and the angular brackets 
denote the ensemble average. The two-point correlation function for the field can be used to analyze, for example, 
self-trapping of light and incoherent soliton propagation in optical structures and nonlinear fibers |l5l Il6l Il7| . Here 
we will use the Fourier transform of the correlation function, as this will give a generalized kinetic description for the 
quasi particles. The Wigner distribution function of the pulse is given by |18|, Il9| 



F(z,t,k,uj) 
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with the inverse F(z+, z_, t+, t_) = J dkduie UT ^F(z, t, k, iv) such that 



I(z,t) 



dk dcu F(z, t, k, u>) 



(7) 



Thus, from Eq. (JJJ the evolution equation for the Wigner function © corresponding to the envelope field E is (see 
also Ref. 
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where we have made a Wigner transformation over the time domain. Here the arrows denote the direction of operation, 
and the operator functions are defined in terms of their respective Taylor expansion. The system of equations Q and 
© determines the evolution of short partially incoherent pulses in nonlinear meta-materials. 

In the low- frequency limit, we retain only the lowest order terms in the sin and cos operators in Eq. JSJ, to obtain 
the Vlasov-like equation 



aud t F - d z F + Kl d t [(l - K 2 I)I]d u F + a [\d t (IF) - lu^I^F] = 0, 



(9) 
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In order to analyse the modulational instability and the effects of the terms due to nonzero a and K2, we write 
F(z, t, uj) = F (ui) +Fx(u>) exp(iKz — iClt) +c.c., where c.c. denotes the complex conjugate, and |.Fi| -C Fq. Expanding 
Eq. JHJ in terms of this ansatz, and using Eq. JjJ, we obtain 



1 = 



{ - Kalo) - a{uj + n/2)] F (u> - 0/2) 
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where Iq = j dujFo(oj). Equation (|10fl represents the NDR for a short electromagnetic pulse, where the pulse may 
have spectral broadening and partial incoherence. 

In the case of a mono-energetic pulse, we have Fq(uj) — IqS(lu — n ), where n corresponds to a frequency shift of 
the background plane wave solution, and the NDR (|1()|> gives 



k = -{2ai + a n )n±n 
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In practice however, the wave envelope will always suffer from perturbations due to various noise sources, e.g. 
thermal spread in the material. A noisy environment may cause the pulse field to attain a random component in its 
phase. Thus, if the phase <p(x) of the electric field varies stochastically, such that the ensemble average of the phase 
satisfies 211 



(exp[— iip(t + t/2)} exp[iip(t — r/2)]) = exp(— firM), 
the background Wigner distribution is given by the Lorentzian spectrum 
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where S7t corresponds to the width of the spectrum. Then, the NDR l|10[l takes the form 



2a [k + ai a n + f n(n - in T )] + aKl (i - n 2 i )n 
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(15) 



This solution generalizes (|ll|l to the case of a random phase background envelope field. The expression Ijl5|) clearly 
shows that the width gives a nontrivial contribution to the NDR. We note that when a = 0, we may define the growth 
rate k according to K = — a^n — iT, and the width then gives rise to a Landau like damping due to i|15|) . 

We may normalize the dispersion relation l|15fl . assuming the GVD to be negative, i.e. a > 0, according to Iq — > kiIq, 
f2x — > s/oinT, n — > ^fan, a — > a/{Ki\fa), and k = k%/ki — \/2nj3. The dispersion relation (|15|l then reads 

k = - (2o-/ + ft - in T ) n ± n [(« + a 2 )i 2 -i + \n 2 + a(n Q - in T )i ] 1/2 ■ (16) 

Letting K = KeK — iT, where Re denotes the real part, we can solve Eq. I|16|) for the growth rate T. Thus, the general 
perturbation exhibits oscillatory growth or damping. The effect of a finite f2-r is in general to reduce the growth rate, 
but it may in certain cases broaden the spectral instability region. We note that a and k have both positive and 
negative values, depending on whether the material is a normal bulk medium, or if it allows for a negative index of 
refraction. For the case of a normal bulk material, we have n > 0, and the higher order nonlinear correction will thus 
reduce the growth rate ^4|- However, for metamaterials we have n < 0, altering the dynamics of Eq. (JTJ in nontrivial 
ways. Thus, for metamaterials, the modulational instability growth rate will differ in measurable ways from the case 
of normal bulk media. We also note that for the modulational instability growth rate, the sign of a is irrelevant. Since 
the full solution is symmetric with respect to an, requiring n > makes the choice of sign for a unimportant in this 
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case. In Fig. ^ we have plotted the growth rate for both positive and negative values of k, so that the difference can 
clearly be seen. Thus, it is reasonable to conclude that soliton formation and filamcntation of ultra-short pulses can 
be much more effective in metamaterials as compared to normal bulk media. 

In fact, using the techniques of previous work (e.g. p3. 123L 123. 125L l26l I27I l28l one can obtain a novel exact 
analytical solution to Eq. representing bright solitons in the form 

E(t, z) = — ^^^^^^=^^^^^^= exp < id arctanh 
J a + cosh 2 [k(vt — z + z Q )] { 

where (p is a constant phase and zq is a spatial displacement, and we require a > 0. Furthermore 
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fl = , a = , Q 5 ' - -, and E% = : , (18) 

2avk^a(l + a) aV 2a 2 k 2 v 3 2' A v 1 

where A = — aKiv[8K2Ct 2 k 2 v 3 — 3(ki<xv + 2a)] + 3cr 2 (l + a 2 k 2 v A ). For a^Owe obtain the solitary waves found in 
Ref. I30L if the material parameters satisfy a 2 = 8aKiK2/3, while a^> 1 corresponds to a flat-top wave form (see Fig. 
EJ- We also not that a nonzero value of a requires a travelling soliton solution (see Fig. [2J ■ 

The width of the pulse is characterized by 1/k. It is however limited by the slowly varying envelope approximation 
used in deducing the basic Eq. (1), and limits how short the pulse can be |3l). This means that the parameter k has 
to be much smaller than the pump wave number feo- 

The issue of stability of solitary solutions is an interesting but complex topic (see e.g. [32I H3)- Bright soliton 
solutions to the one-dimensional cubic nonlinear Schrodinger equation are known to be unstable to transverse per- 
turbations, but it is also known that competing nonlinearities, occuring in cubic-quintic systems, can stabilize soliton 
solutions 34]. Our equation contains such cubic-quintic nonlinearites, as well as a derivative nonlinearity. It is 
therefore not unreasonable to conjecture that the solution <|17fl is stable. However, in order to strictly answer the 
question of stability of l|17|l much further research is needed. 

To summarize, we have examined the modulational instability and localization of an ultra-short electromagnetic 
pulse that is governed by a modified nonlinear Schrodinger equation, which is also valid for negative refractive index 
media. We have considered both coherent and incoherent pulses and derived nonlinear dispersion relations that depict 
the growth rates of the modulational instabilities. The effect of partial coherence of the electromagnetic pulse is to 
reduce the instability growth rate. Furthermore, we have found that modulationally unstable pulses can appear as 
a new type of light envelope solitons whose profiles are found analytically and analyzed numerically. Localized light 
pulses can propagate without distortion and should be observable in metamaterials which have negative index of 
refraction. 
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FIG. 1: The modulational instability growth rate plotted for Io = 0.5, using a = 0.5, Qt = 0.1, and k = —0.5 (upper curve) and 
ft = 0.5 (lower curve). We have here normalized the growth rate by Kilo and the frequency by v/|/ii|/o/<ai, and correspondingly 
for the parameters in Eq. 115L thus giving dimensionless quantities. The larger growth rate for metamaterials is clearly seen. 
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FIG. 2: The characteristic features of the soliton solution Q17[l . The real part of the solution 1)17^ is plotted in panel (a), while 
the intensity / = \E\ 2 is plotted in panel (b), both for two different values of z. We have normalized the solution so that 
E — > ^JaE/Eo, while time is normalized by ^/\h\\Eq / a. 



